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Abstract — A new model for a cluster of hybrid sensors network 
with multi sub-clusters is proposed. The model is in particular 
relevant to the early warning system in a large scale monitoring 
system in, for example, a nuclear power plant. It mainly addresses 
to a safety critical system which requires real-time processes with 
high accuracy. The mathematical model is based on the extended 
conventional search algorithm with certain interactions among the 
nearest neighborhood of sensors. It is argued that the model could 
realize a highly accurate decision support system with less number 
of parameters. A case of one dimensional interaction function is 
discussed, and a simple algorithm for the model is also given. 

I. Introduction 

A safety critical system in, for example, nuclear power plants 
involves complicated and advance safety system for ensuring 
its running process absolutely safe. The system consists of 
high number and various types of sensors (T|, 0. All of them 
generate a huge amount of data at a real-time basis which must 
be processed properly throughout its whole life time. Many 
heuristic approaches based on the artificial intelligence (AI) 
such as neural network are currently available and have been 
studied intensively [3 |-[6|. However, the Al-based approaches 
have a fundamental problem due to its statistical algorithm 
which could lead to disaster in the real applications of safety 
critical system. 

It is clear that in such critical systems, no fault tolerance 
is the most important principle. Therefore, putting the safety 
as the priority, one should implement the exhaustive algorithm 
spanning over all possibilities rather than using the AIs. This 
turns into the exhaustive search problem which unfortunately 
lacks of inefficiency due to the requirement of huge computing 
power. Some technical approaches have been introduced to 
overcome this problem. Most of them deploy the parallel 



algorithm J7), (8] together with graphical or combinatorial 
representation to improve both resources and running time [9 1- 

m. 

Previously, the application of exhaustive methods into search 
problem was not feasible for large number of sensors. For- 
tunately, the affordable parallel environment using graphic 
processor unit (GPU) is available in recent days. The use of 
GPU is getting popular, especially, after the introduction of the 
NVIDIA Compute Unied Device Architecture (CUDA) through 
a C-based API [ 13 1. This enables an easy way to take advantage 
of the high performance of GPUs for parallel computing. De- 
ploying the GPU-based distributed computing would reduce the 
execution time causing less responsive system in the previous 
days, while it also realizes lower power and space consumption 
than CPU [14|. This motivates us to reconsider the feasibility 
of exhaustive search for hybrid sensors network. 

In this paper a new exhaustive search based model is pro- 
posed. The model is mainly intended to realize an exhaustive 
decision support system (DSS) consisting of various and huge 
number of sensors. However, the paper is focused only on 
introducing the model. The discussion of parallelization and 
detail analysis will be published elsewhere. 

The paper is organized as follows. First in Sec. [II] the 
model is introduced, and it is followed by the description 
of mathematical formulation in Sec. [Ill] In Sec. [IV] a simple 
algorithm to execute the model is given. Finally, the paper is 
ended with a short summary and discussion. 

II. The model 

Before moving on constructing the model and its mathe- 
matical representation, let us discuss the basic constraints and 



circumstances in the expected applications. Putting in mind that 
the model is developed under the following considerations : 

1) No failure decision is allowed. 

By definition, there is no room for even a small mistake 
generated by the DSS. This actually discourages the 
deployment of any Al-based method from first principle. 

2) Fast enough 'real-time' process 

Fast execution time of the whole process is crucial to 
increase the safety. However, fast execution time at the 
order of minutes is in practical more than enough. This 
moderate requirement encourages the implementation of 
exhaustive methods supported by GPU powered compu- 
tation. 

3) Huge number of hybrid sensors 

The system is consisting of huge number (at the order of 
hundreds or thousands) sensors with different character- 
istics, in particular types and scales fl5l . 

4) Certain relationship across the sensors 

Each sensor has relations with another ones in certain way 
at some degrees. The relationships among the sensors are 
thereafter called as interaction. 

5) Sensor network with multi clusters 

The sensor network is divided into several clusters which 
typically represents the geographical locations with dif- 
ferent degrees of interaction. This realizes a situation of, 
for instance, a nuclear power plant which is equipped 
with many sensors in each building. Consequently the 
sensors in each cluster have stronger interactions among 
themselves than with another ones belong to another 
clusters. So, the model should be able to describe the 
independent interactions among the sensors in a cluster, 
and also the interactions among different clusters as well. 

6) Dynamic behavior 

The values of each sensor are by nature changing from 
time to time. However, the data acquisition is performed 
periodically, for instance every few minutes according 
to the above second point. In a nuclear reactor facility 
this could happen due to human errors, common system 
failures and even seismic activities Ifl6ll . 

Having the above requirements in mind, obviously one will 
arrive at the problem of unlimited decision trees. In order to 
reduce the tree significantly without raising the risk of failures, 
let us assume the nearest neighborhood approximation (NNA). 
Under this approximation, only the interactions with the nearest 
sensors are taken into account. 

In the present paper, let us consider the simplest case of one 
dimensional (1-D) relationship. This means all sensors are put 
on a virtual line which allows only interactions with the nearest 
right and left neighboring sensors for each sensor. This actually 
reproduces the known tree analysis commonly implemented in 
the analysis of fault Q6), Q7), elements interaction Ifl8l-ll20l 
and even in optimizing system ETL 

The algorithm can be well illustrated in a tree-like diagram 
using the evaluation scheme under the NNA scheme as shown 
in Fig.Q] In the figure, two adjacent sensors are first evaluated 




Fig. 1: The illustration of 1-D relationship of n sensors eval- 
uated till the m th level, where S denotes sensor and E is the 
evaluation result. 



and then the result are subsequently evaluated with another 
adjacent sensor. The evaluation scheme depicted in Fig. Q] can 
be exhaustively changed according to the acquired values from 
the responsible sensors. It should be noted that the evaluation 
scheme is not necessarily binary, but it could be anything else 
like fuzzy and so on. 

On the other, due to point 4 one should consider the modified 
tree analysis, that is both edges should interact each other too 
and forms a circle line of sensors. Moreover, each sensor on 
the circle line should be put carefully. Because, according to 
point 5 and the NNA, the relative location of sensors on the 
circle line represents their degree of relationship or relevancy 
between one and another. The stronger relationship between 
two sensors, both should be put closer each other. This type of 
circular model is depicted in Fig. [2] There are two examples 
of evaluation results in the figure, the blue and brown ones 
corresponding to the evaluation at different time. The innermost 
circles represent the chain of sensors, and the subsequent outer 
circles describe the evaluation results at certain levels. 

As required in point 5, the sensor network should also 
be divided into several clusters based on either its genuine 
characteristic or critical levels. Each cluster can be treated 
separately as an independent sensor network as Fig. [2] The 
model with several clusters is illustrated in Fig. [3] where each 
cluster is separated by the blue dashed lines. Two evaluation 




Fig. 2: The circular representation for 1-D relationship of n 
sensors evaluated for two schemes (blue and brown lines) till 
the m th level, where S denotes sensor and E is the evaluation 
result. 



results are shown in the figure as before, the blue and brown 
ones corresponding to the evaluation at different time. 
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Fig. 3: The circular representation for 1-D relationship of n 
sensors belong to separated clusters, evaluated for two schemes 
(blue and brown lines) till the m th level, where S denotes 
sensor and E is the evaluation result. 



TABLE I: The values of weight parameter for n sensors. 
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Now, one is ready to formulate the model in a mathematical 
representation. 

III. Mathematical representation 

Based on the previous discussion, one should first model the 
interaction between the neighboring sensors. This should be a 
function which determines the value of Eij representing the 
interacting result between two adjacent sensors at the point i 
and j, 

1 ' ! ' 1 (1) 



Eij — &ij 



k+1 U+l 



Here, k — 1, 2, • • • , m denotes the evaluation level of Xi with 
i = 1, 2, • • • , n and m = n — 1. While xi is the normalized 
value acquired by sensor for U = 1, or the value of previous 
evaluation result for li > 1. The coupling constant ay reflects 
the strength of interaction between both adjacent sensors, Xi 
and Xj respectively. It is defined in a way such that the value 
of Eij is getting smaller for higher evaluation level, that is, 
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(2) 



for li = 1 in a system with n sensors. 

Moreover, Eqs. ([TJ and © can be extended for > 1 as 
follow, 

Xij X]^ 
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with, 
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(3) 



(4) 



and x^ = E^ respectively. The situation is well illustrated in 
Figs. |2] and [3] Further generalization of Eqs. (fJJ and (0 up to 
certain evaluation level can be performed in a straightforward 
way. 

There is actually an important reason for choosing the 
definition in Eqs. (f2]i and (0J. In order to fulfill the requirements 
in point 3 and 4 in Sec. [TTJ it is plausible to normalize all 
scales in a uniform unit scale. In the present case all values 
are normalized to be in the range between and 1. This 
normalization demands all acquired values from the sensors 
should also be normalized accordingly through a relation, 







(5) 



where X{ and x\ are the normalized and originally acquired 
values of i th sensor for U = 1. The normalization factor / is, 

1 



(6) 



with £ m ax/min denotes the maximum or minimum value of each 
sensor. This kind of normalization enable us to treat all sensors 
in the same manner regardless with its types and unit scales. 

From Eqs. (0 and (0J, it is obvious that the coupling a 
is uniquely characterizing the present model. It ensures the 
evaluation value of E at the final (n — l) th level is always 
divergent, that is between and 1. This is in contradiction 
with any conventional tree analysis which associates the largest 
evaluation value at the final level as the final solution. By the 
way, from the definition 1/n < ay < 1 as shown in Tab. U 
which forms a symmetric matrix with unit diagonal elements. 

Furthermore, one should take a threshold value Eth as a 
standard value whether the evaluation value at certain level is 
allowed to proceed further or not. For the sake of simplicity, 
this value is fixed and valid for all levels and sensors. This 
represents the critical value of safety. Following the above 
normalization it is again constrained, 



Algorithm QJ3] require the array of sensor's indexes and 
values. The value described in a sensor's index array absolutely 
points to its value in the array of sensor's value. Algorithm [T] 
determines the positioning of sensors, while Algorithm [2] and 
[3] evaluate their interactions. 

Algorithm 1 Positioning sensors 

Require: root {index of sensor being a root of a tree or a sub 
tree} 

Require: n {number of sensors involved} 
Require: index {array of sensor's index} 
1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



< E th < 1 



(7) 



Only if the evaluation value exceeds this threshold, i.e. E > 
E t h, the tree should be analyzed further. Otherwise it ends 
forever. According to the initial value of each sensor at certain 
time, some evaluation values at final level may survive or 
not. The surviving value triggers the warning alarm which is 
indicating some anomalies detected by any sensors. Of course, 
the determination of appropriate E t h requires preliminary ex- 
periments based on the available standards and regulations. 

The above procedure is carried out each time following the 
periodic data acquisition by all sensors. Finally, all the tools 
have been established and we are ready for applying the above 
rules. 

IV. The algorithm 

In this section, let us provide a simple algorithm to realize 
the previously discussed model. 

Because of the exhaustive method deployed in the model, 
it requires all variables are treated as a circle. The array of 
variables should be a traversing pointer indicating the first 
element as a starting point. Further, the evaluation is performed 
from the starting point, proceeds to the subsequent element of 
variables array in an increasing mode till reaching the last one, 
that is back to the first element. Each sensor has a chance to 
become a root of a tree and also a leaf. 

While forming a circle model, the tree within the model is 
evaluated recursively. A set of simple algorithms is presented 
here. It consists of two main parts for the pra-evaluation and 
the main evaluation. Each sensor is labeled with an integer 
ranging from to n — 1 for n number of sensors. The algorithm 
starts with positioning the sensors by shifting it one by one to 
generate considerable combinations. Each time a tree of sensors 
is formed, it is evaluated as a scheme depicted in Fig. Q] 



if n = 2 or root — n — 2 then 
evaluate their interaction 

exchange sequence of the two or the last two sensors 
evaluate their interaction 
else 

for i = — > n do 

re-positioning sensorsfmdea;, root + 1) 
shifting the sequence from the root position 
end for 
end if 



Algorithm 2 Calculating interaction's weight 

Require: root {index of sensor being a root of a tree or a sub 
tree} 

Require: n {number of sensors involved} 
Require: index {array of sensor's index} 
Require: w = {weight of sensor's interaction initialized} 
Require: 6 = {counting number of sensors combination 
currently involved in interaction} 
for i = 1 — >■ root do 



for j = i 
w = w 
& = 64 
end for 
end for 
return^ 



(1- 



root do 

abs(index[i\- 



■index[j\) 



Algorithm 3 Evaluate the interaction 



Require: x\, X2 {two sensor value involved, for a more deeper 
tree, xl can be an evaluation value from the previous level 
of depth} 

Require: h,fa {level of interaction} 

Require: a {weight of interaction} 

!■ ^-"((irrr) + (i^T)) 

In positioning the sensors as determined by Algorithm Q] 
for the case of two sensors the algorithm only exchanges their 
sequences. After the evaluation of first sequence, the algorithm 
exchanges the sequence and reevaluate the new one. This 
procedure is specified in the second to fourth lines. For the case 
of more than two sensors, the algorithm recursively traverses 



the sequence till it reaches the condition where only two sensors 
are left. In this case, the algorithm runs in the same manner as 
if there are only two sensors involved. Each time the algorithm 
is traversing deeper, the root is increased indicating the depth of 
tree under evaluation. If a leaf of tree is reached, the algorithm 
returns back to the parent leaf, exchanges to the next root, 
traverses deeper till it reaches the leaf. After all paths to each 
leaf have been reached by one root, the sequence exchanges to 
another sensor as the new root and so forth. 

Furthermore, the weight of interaction is calculated using 
Eq. d2J and given in Tab. U In Algorithm Q]roo£ determines the 
root of a new sub tree, while in Algorithm [2] root is intended to 
determine the number of sensors which are currently involved 
in the interaction as illustrated in Fig. __ b in Algorithm __ is 
intended to count the number of combination of two sensors 
among the whole sensors as defined in Eq. __ that is it would 
be as many as elements in Tab. |U For example, the interaction 
of three sensors contains three dual-sensor interactions, while 
the interaction of four sensors contains six and so on. 

Finally, Algorithm [3] is intended to calculate the evaluation 
value, E. It requires the weight of each interaction from 
Algorithm __ using Eq. __) and Tab. _] The total number of 
evaluation values is equal to the number of sensors being 
involved. 

V. Summary 

A new model based on the exhaustive search method for 
hybrid sensor network has been proposed. The model treats 
all sensors in the same manner by introducing normalization 
procedure for all sensors and parameters. It is shown that the 
model is able to describe the whole evaluation processes using 
few parameters, that is the coupling constant a for each pair of 
sensors determined uniformly and the universal threshold value 

In the present paper, the study is focused on the case of 
sensor network with 1-D relationship. A simple algorithm 
for such cases has also been given and briefly discussed. 
Through the discussion, it is argued that the model could 
realize a feasible early warning system for any safety critical 
facilities involving various sensors using exhaustive method to 
prevent unnecessary failures due to statistical approaches in, 
for example any Al-based methods. On the other hand, the 
method requires much less computing power since it has only 
a complexity of 0(n\). 

In principle the method can be extended to incorporate more 
complicated relationship among the sensors by considering 
higher dimensional relationship. Nevertheless, some studies on 
distributing the computation load to improve the processing 
speed should also been done carefully. All of these issues are 
in progress and will be published elsewhere. 
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